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Abstract 

Proton-proton collision data recorded in 2011 and 2012 by the LHCb exper¬ 
iment, corresponding to an integrated luminosity of 3.0 fin 1 , are analysed to 
search for the charmless B° —» decay. More than 600 B° —>• (7r + 7r _ )(7r + 7r _ ) 

signal decays are selected and used to perform an amplitude analysis, under the 
assumption of no CP violation in the decay, from which the B° —> p°p° decay 
is observed for the first time with 7.1 standard deviations significance. The 
fraction of B° —> p°p° decays yielding a longitudinally polarised final state is 
measured to be /l = 0.745]lo'^(stat) ± 0.034(syst). The B° — > p°p° branch¬ 
ing fraction, using the B° —> 4>K*( 892)° decay as reference, is also reported as 
B{B° -> p°p°) = (0.94 ± 0.17(stat) ± 0.09(syst) ± 0.06(BF)) x 10“ 6 . 
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1 Introduction 


The study of B meson decays to pp final states provides the most powerful constraint to 
date for the Cabibbo-Kobayashi-Maskawa (CKM) angle a = arg [{V t .dV t l)/(V u dV* b )\ jT]{3]. 
Most of the physics information is provided by the decay B° —y p + p~ as measured at the 
e + e~ colliders at the T(4S) resonance 00 Q for which the dominant decay amplitude, 
involving the emission of a W boson only (tree), exhibits a phase difference that can be 
interpreted as the sum of the CKM angles (3 + 7 = n — a in the Standard Model. The 
subleading amplitude associated with the exchange of a W boson and a quark (penguin) 
must be determined in order to interpret the electroweak phase difference in terms of the 
angle a. This is realised by means of an isospin analysis involving the companion modes 
B + —>• p + p° |6,7| and L>°—>• p°p° [8,91 j^] In particular, the smallness of the amplitude of 
the latter leads to a better constraint on a. 

The BaBar and Belle experiments reported evidence for the B° —> p°p° decay |8][9| 
with an average branching fraction of B(B° —> p°p°) = (0.97 ± 0.24) x 10 -6 ||8,|9|. Despite 
small observed signal yields, each experiment measured the fraction /l of decays yielding a 
longitudinally polarised final state through an angular analysis. The Belle collaboration did 
not find evidence for polarisation, /l = 0.211q 26 |9]> while the BaBar experiment measured 
a mostly longitudinally polarised decay, /l = 0.75Iq!i 5 j&j- These results differ at the level 
of 2.0 standard deviations. The large LHCb data set may shed light on this discrepancy. 
In addition, LHCb may confirm the hint of B° —» p°/ 0 (980) decays reported by Belle |9|. 
Measurements of the B° —>■ p°p° branching fraction and longitudinal polarisation fraction 
at LHCb can be used as inputs in the determination of a 12,3]. 

This work focuses on the search and study of the B° —> (7r + 7r _ )(7r + 7r _ ) decay in which 
the two (7r + 7r _ ) pairs are selected in the low invariant mass range (< 1100 MeV/c 1 2 ). The 
B° —)■ p°p° is expected to dominate the (7r + 7r“) mass spectrum. The (7r + 7r“) combinations 
can actually emerge from S-wave non-resonant and resonant contributions or other P- or 
D-wave resonances interfering with the signal. Hence, the determination of the B° —> p°p° 
yields requires a two-body mass and angular analysis, from which the fraction of the 
longitudinally polarised final state can be measured. 

The branching fraction is measured relative to the B° —>■ 892)° mode. The 

B° —> 4>K*( 892)° decay, which results in four light mesons in the final state, is similar to 
the signal, thus allowing for a cancellation of the uncertainties in the ratio of selection 
efficiencies. 


2 Data sets and selection requirements 

The analysed data correspond to an integrated luminosity of l.Ofb -1 and 2.0fb -1 from 
pp collisions recorded at a centre-of-mass energy of 7TeV, collected in 2011, and 8TeV, 
collected in 2012, by the LHCb experiment at CERN. 

1 Charge conjugation is implicit throughout the text unless otherwise stated. 

2 p° stands for p°(770) throughout the text. 
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The LHCb detector @0 is a single-arm forward spectrometer covering the 
pseudorapidity range 2 < 77 < 5, designed for the study of particles containing 6 or 
c quarks. It includes a high-precision tracking system consisting of a silicon-strip vertex de¬ 
tector surrounding the pp interaction region 12 , a large-area silicon-strip detector located 


upstream of a dipole magnet with a bending power of about 4 Tm, and three stations of 
silicon-strip detectors and straw drift tubes 13 placed downstream of the magnet. The 
tracking system provides a measurement of momentum, p, of charged particles with a 
relative uncertainty that varies from 0.5% at low momentum to 1.0% at 200GeV/c. The 
minimum distance of a track to a primary vertex, the impact parameter, is measured with 
a resolution of (15 + 29/px) pm, where p t is the component of the momentum transverse 
to the beam, in GeV/c. Different types of charged hadrons are distinguished using in¬ 
formation from two ring-imaging Cherenkov (RICH) detectors 114]. Photons, electrons 
and hadrons are identified by a calorimeter system consisting of scintillating-pad and 
preshower detectors, an electromagnetic calorimeter and a hadronic calorimeter. Muons 
are identified by a system composed of alternating layers of iron and multiwire proportional 
chambers 15 . The online event selection is performed by a trigger 16 , which consists of 


a hardware stage, based on information from the calorimeter and muon systems, followed 
by a software stage, which applies a full event reconstruction. 

In this analysis two categories of events that pass the hardware trigger stage are 
considered: those where the trigger decision is satisfied by the signal 6 -hadron decay 
products (TOS) and those where only the other activity in the event determines the 
trigger decision (TIS). The software trigger requires a two-, three- or four-track secondary 
vertex with large transverse momenta of charged particles and a significant displacement 
from the primary pp interaction vertices (PVs). At least one charged particle should 
have Pt > 1.7 GeV/c and is required to be inconsistent with originating from any primary 
interaction. A multivariate algorithm 17 is used for the identification of secondary vertices 
consistent with the decay of a 6 hadron. 

Further selection criteria are applied offline to reduce the number of background events 
with respect to the signal. The ( 7 r + 7 r _ ) candidates must have transverse momentum larger 
than 600MeV/c, with at least one charged decay product with px > lOOOMeV/c. The two 
( 7 T+ 7 T - ) pairs are then combined to form a B° candidate with a good vertex quality and 
transverse momentum larger than 2500 MeV/c. The invariant mass of each pair of opposite- 
charge pions forming the B° candidate is required to be in the range 300-1100 MeV/c 2 . 
The identification of the final-state particles (PID) is performed with dedicated neural- 
networks-based discriminating variables that combine information from the RICH detectors 
and other properties of the event 14 . The combinatorial background is further suppressed 


with multivariate discriminators based on a boosted decision tree algorithm (BDT) [18.19 
The BDT is trained with simulated B° —» p°p° (where p° —>■ 7 r + 7 r _ ) events as signal 
sample and candidates reconstructed with four-body mass in excess of 5420 MeV/c 2 as 
background sample. The discriminating variables are based on the kinematics of the B 
decay candidate (B px and the minimum px of the two p° candidates) and on geometrical 
vertex measurements (quality of the B candidate vertex, impact parameter significances 
of the daughters, B flight distance significance and B pointing to the primary vertex). 
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The optimal thresholds for the BDT and PID discriminating variables are determined 
simultaneously by means of a frequentist estimator for which no hypothesis on the signal 

The 77° meson candidates are accepted in the mass range 5050- 
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yield is assumed 
5500 MeV/c 2 . 

The normalisation mode 77°—)- <pK*( 892)° is selected with similar criteria, requiring 
in addition that the invariant mass of the ( K + ti ~) candidate is found in a range of 
±150 MeV/c 2 around the known value of the K* (892)° meson mass 21 and the invariant 
mass of the (. K + K ~) pair is in a range of ±15 MeV/c 2 centred at the known value of the 
0 meson mass |21|. A sample enriched in 77° —» (/i + 7r~)(7T + 7r _ ) events is selected using 
the same ranges in (7r + 7r“) and (K + 7i~) masses to estimate the background with one 
misidentified kaon. 

The presence of (7r + 7r“) pairs originating from ,7/0, Xco and Xci charmonia decays is 
vetoed by requiring the invariant masses M of all possible (7r + 7r“) pairs to be | M — M 0 \ > 
30 MeV/c 2 , where M 0 stands for the corresponding known values of the J/0, Xco and Xci 
Similarly, the decays D° —> K~n + and D° —> 7r + 7r - are vetoed by 


meson masses 
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requiring the corresponding invariant masses to differ by 25 MeV/c 2 or more from the 
known 77° meson mass [2l|. To reduce contamination from other charm backgrounds and 


from the 77°—a0(—>■ p°7r + )7r _ decay, the invariant mass of any three-body combination 
in the event is required to be larger than 2100 MeV/c 2 . 

Simulated 77°-)- p°p° and 77°—)- 4>K*( 892)° decays are also used for determining the 


relative reconstruction efficiencies. The pp collisions are generated using PYTHIA 
with a specific LHCb configuration 1231. Decays of hadronic particles are described by 
EvtGen 24 . The interaction of the generated particles with the detector and its response 


are implemented using the Geant4 toolkit 25 as described in Ref. 26 


3 Four-body mass fit 

The four-body mass spectrum M(7r + 7r _ )(7r + 7r _ ) is fit with an unbinned extended likeli¬ 
hood. The fit is performed simultaneously for the two data taking periods together with 
the normalisation channel M(K + K~)(K + n~) and PID misidentification control channel 
mass spectra. The four-body invariant mass models account for B° 
and possible 77/ signals, combinatorial backgrounds, signal cross-feeds and background 
contributions arising from partially reconstructed 6-hadron decays in which one or more 
particles are not reconstructed. 

The 77° and 77/ meson shapes are modelled with a modified Crystal Ball distribution 
A second power-law tail is added on the high-mass side of the signal shape to account 
for imperfections of the tracking system. The model parameters are determined from 
a simultaneous fit of simulated signal events that fulfill the trigger, reconstruction and 
selection chain, for each data taking period. The values of the tail parameters are identical 
for the 77° and 77/ mesons. Their mass difference is constrained to the value from Ref. 

The mean and width of the modified Crystal Ball function are free parameters of the fit 
to the data. 
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Figure 1: Reconstructed invariant mass spectrum of (left)(7r + 7r~)(7r + 7r _ ) and (right)(/v + Jv _ )(A' + 7t _ ). 
The data are represented by the black dots. The fit is represented by the solid blue line, the B° signal by 
the solid red line and the B ° by the solid green line. The combinatorial background is represented by the 
pink dotted line, the partially reconstructed background by the cyan dotted line and the cross-feed by the 
dark blue dashed line. 


The combinatorial background in each four-body spectrum is described by an exponen¬ 
tial function where the slope is allowed to vary in the fit. 

The misidentihcation of one or more final-state hadrons may result in a fully recon¬ 
structed background contribution to the corresponding signal spectrum, denoted signal 
cross-feed. The magnitude of the branching fractions of the signal and control modes 
as well as the two-body mass selection criteria make these signal cross-feeds negligible, 
with one exception: the misidentihcation of the kaon of the decay B° —>■ (Ji + 7r _ )(7r + 7r _ ) 
as a pion yields a significant contribution in the A4(7t + 71~)(tt + 71~) mass spectrum. The 
mass shape of B° —> (K + 7t~)(7t + 7t~) decays reconstructed as B° — > (7r + 7T~)(7T + 7i~) is 
modelled by a Crystal Ball function, whose parameters are determined from simulated 
events. The yield of this signal cross-feed is allowed to vary in the fit. The measurement 
of the actual number of reconstructed B° — > (K events multiplied by the 
data-driven estimate of the misidentihcation efficiency is consistent with the measured 
yield. 

The partially reconstructed background is modelled by an ARGUS function 
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convolved with a Gaussian function accounting for resolution effects. Various mass shape 
parameterisations are examined. The best ht is obtained when the endpoint of the ARGUS 
function is hxed to the value expected when one pion is not attributed to the decay. The 
other shape parameters of the ARGUS function are free parameters of the ht, common to 
the two data taking periods. The floating width parameter of the signal mass shape is 
constrained to be equal to the width of the Gaussian function used in the convolution. 

Figure [l] displays the M(7r + 7r - )(7r + 7r - ) and M(K + K~)(K + n~) spectra with the ht 
results overlaid. The signal event yields are shown in Table [lj Aside from the prominent 
signal of the B° —> () decays, the decay mode B° s —> (7r + 7r - )(7r + 7r - ) is 
observed with a statistical significance of more than 10 standard deviations. The statistical 
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Table 1: Yields from the simultaneous fit for the 2011 and 2012 data sets. The first and second 
uncertainties are the statistical and systematic contributions, respectively. 


Decay mode 

Signal yields 2011 

Signal yields 2012 


+ 

A 

+ 

i 

185 ± 15 ± 4 

449 ± 24 ± 

7 

B° 

+ 

A 

+ 

1610 ± 42 ± 5 

3478 ± 62 ± 10 

B°-> 

(K + K~ ) (K + tt~ ) 

1513 ±40 ±8 

3602 ± 62 ± 10 


+ 

i 

A 

+ 

30 ± 7 ± 1 

71 ± 11 ± 

1 

B ° s -> 

+ 

40 ± 10 ± 3 

96 ± 14 ± 

6 

B ° s -> 

+ 

+ 

42 ± 10 ± 3 

66 ± 13 ± 
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significance is evaluated by taking the ratio of the likelihood of the nominal fit and of the 
fit with the signal yield fixed to zero. 

A systematic uncertainty due to the fit model is associated to the measured yields. 
The dominant uncertainties arise from the knowledge of the signal and signal cross¬ 
feed shape parameters determined from simulated events. Several pseudoexperiments 
are generated while varying the shape parameters within their uncertainties, and the 
systematic uncertainties on the yields are estimated from the differences in results with 
respect to the nominal fit. 


4 Amplitude analysis 


An amplitude analysis is used to determine the vector-vector (VV) contribution B° —> p u /F 
by using two-body mass spectra and angular variables. The four-body mass spectrum is 


first analysed with the sPlot technique 29 to subtract statistically the background under 
the B° —>■ (7r + 7r _ )(7T^ 


7T 


signal. 


For the two-body mass spectra, contributions from resonant and non-resonant scalar 
(S'), resonant vector ( V ) and tensor (T) components are considered in the amplitude 
fit model through complex mass propagators, M(?7p), where the label i — 1,2 are the 
first and second pion pairs, which are assigned randomly in every decay since they are 
indistinguishable. The P-wave lineshape model comprises the p° meson, described using 
the Gounaris-Sakurai parameterisation 30 , and the c o meson, parameterised with 

a relativistic spin-1 Breit-Wigner M w (m,j). The D-wave lineshape Mf 0 (rrii) accounts for 
the /"2 (1270), modelled with a relativistic spin-2 Breit-Wigner. The S-wave model includes 
the /o(980) propagator Mj( 980 )(mj), described using a Flatte parameterisation |31,32], 
and a low-mass component. The latter includes the broad low-mass resonance /o(500) 
and a non-resonant contributions, which are jointly modelled in the framework of the 
A"-matrix formalism 33 and referred as M( 7r7r ) 0 (mj). Following the A'-matrix formalism, 
the amplitude for the low-mass 7r + 7r _ S-wave can be written as 


A(m) oc 


K 


1 - ipK 


( 1 ) 
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with 


K = K + K 

11 — -*^res i non—res 


m 0 r(m) 


(mg — m 2 )p{m ) 


+ K, 


p(m) = 2 


g(m) 


m 


( 2 ) 

( 3 ) 


where k is measured to be —0.07 ± 0.24 from a fit to the inclusive 7 r + 7 r - mass distribution 
and mo and T are the nominal mass and mass-dependent width of the / 0 (500), as determined 
in Ref. 34 . The functions p(m) and q{m ), defined in Ref. 33 , are the phase space factor 


and the relative momentum of a pion in the p° centre-of-mass system. By convention, the 
phase of the iV /( 7r7r ) 0 (rnf) mass propagator is set to zero at the p° nominal mass. 

The signal sample is described by considering the dominant amplitudes of the signal 
decay. The B —» VV component contains the B —> p°p° and B° —>• p°u) amplitudes. 
The B —* VS component accounts for B° —* p°(tt + tt~) 0 and B° —> p°/ 0 (980) amplitudes 
and the B —> VT contribution is limited to the purely longitudinal amplitude of the 
B° —> p°f 2(1270) transition. Because of the broad natural width of the af particle, a small 
contamination from the decays B° —> afi r T remains in the sample. This contribution with 


,± 


ap —> p 0 7 r ± in S-wave is considered along with its interference with the other amplitudes. 
This is done by introducing the CP-even eigenstate from the linear combination of 
individual amplitudes of the decays B° 


a{ 7 T 


and P u 


a ± 7i n 


as dehned in Ref. 35 


The contribution of the decays B° — > ojoj, B° —> / 0 (980)/o (980), B° —> ojS, B° —>■ ujT, 
B° —> / 2 (1270)S', B° —> / 2 (1270)/2(1270) and B° —> (p°/ 2 ( 1270))|| i j_ are assumed to be 
negligible, where the || and _L subindices indicate the parallel and perpendicular amplitudes 
of the decay. The choice of the baseline model was made prior to the measurement of 
the physical parameters of interest after comparing a set of alternative parameterisations 
according to a dissimilarity statistical test [36]. 

The differential decay rate for B° —y ( 7 r + 7 r - )( 7 r + 7 r - ) decays at the B° production time 
t — 0 is given by 


d 5 T 

d cos 9 1 d cos 9-2 d<p d m\ dm 2 


oc $ 4 (mi,m 2 ) 


^ A fi ( ? «i pm 2y 9 1 ,9 2 ,p) , 

i=l 


( 4 ) 


where the variables 9 1 , 9 2 and <p are the helicity angles, described in Fig. [2j and $4 is 
the four-body phase space factor. The notations of the complex amplitudes, Ai, and the 
expressions of their related angular distributions, fi, are displayed in Table [2] The mass 
propagators included in the fi functions are normalised to unity in the fit range. 

For the CP conjugated mode, B° —$■ the decay rate is obtained under 

the transformation Ai —> ipAi, where r/* is the CP eigenvalue of the CP eigenstate i, shown 
in Table [U 

The untagged time-integrated decay rate of B° and B° to four pions, assuming no CP 
violation, can be written as 


_d 5 (r + r)_ 

d cos 9\ d cos 9 2 dp d m\ dmf 


« 2 J 2 J Ke[4 i ^*/ i /*](2 
j = 1 i<j 


Oij){ 1 • yi’ 
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Table 2 : Amplitudes, Ai, CP eigenvalues, 77 *, and mass-angle distributions, f\ , of the 
B° —> ( TT + TT~)(Tr + 7r~ ) model. The indices ijkl indicate the eight possible combinations of pairs 
of opposite-charge pions. The angles a^i, ffj and &ki are defined in Ref. 
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Ai rji 


A 

1 

A II 

S±pp 

1 


-1 

A% 

1 

A II 

Sl-puj 

1 

A- L 

Sipcj 

-1 

Ap(7rir)o 

-1 

Apf( 980) 

-1 

^( 7r7r )o( 7r7r )o 

1 

A0 pf2 

-1 


1 


_ fi _ 

M p (toi)M p (to 2 ) cos 8 1 cos 8 2 

Mp{mi)Mp{m 2 sin #1 sin 60 cos 
M p (mi)M p (ro 2 ) ^ sin sin 8 2 sin ^ 

^[M p (mi)M w (m 2 ) + Af w (TOi)M p (m 2 )] cos#i cos# 2 
[M p (mi)M w (m 2 ) + M u) (mi)M p (m 2 )] ^ sin sin 8 2 cos tp 
^ [M p (mi) M u (m 2 ) + M w (mi)M p (m 2 )]^ sin^i sin0 2 sin</? 

^ [Af p (mi)M (7rOT ) 0 (m 2 ) cos 6>i + M (7T7r)o {m\)M p (m 2 ) cos 0 2 ] 
^[M p (mi)M /(980 )(m 2 ) cos#! + A7 /(980 )(mi)M p (m 2 ) cos# 2 ] 

_ Tf ( 7 rjr)o (mi)M( ^) 0 (m 2 )| 

[M p (mi)My 2 (m 2 ) cos#i(3cos 2 # 2 — 1) + M/ 2 (mi)M p (ro 2 ) cos# 2 (3cos 2 #1 — 1)] 
Efafci} 71 ^ai {jnjjk)Mp(mij) [cos a fc ; cos + sin sin f3 ik cos 



Figure 2: Helicity angles for the ( 7 r + 7 r )( 7 r + 7 r ) system. 


where 8 t j = 1 when i — j and 8 i3 = 0 otherwise. 

The efficiency of the selection of the final state B° —y ( 7 r + 7 r _ )( 7 r + 7 r _ ) varies as a function 
of the helicity angles and the two-body invariant masses. To take into account variations 
in the efficiencies, four event categories k are defined according to their hardware trigger 
decisions (TIS or TOS) and data taking period (2011 and 2012 ). 

The acceptance is accounted for through the complex integrals 

Uij = J e(6i,9 2 , v?,mi,m 2 )/j/* (2 - <Sjj)$ 4 (mi, m 2 )dcos0i dcos0 2 d</?dmi dm^, (6) 

where fi are the functions given in Table [2] and e the overall efficiency. The integrals are 
computed with simulated events of each of the four considered categories, selected with 
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Table 3: Results of the unbinned maximum likelihood fit to the angular and two-body invariant 
mass distributions. The first uncertainty is statistical, the second systematic. 


Parameter 

Definition 

Fit result 

/l 

|App| 2 /(|A° p | 2 + |App 2 + l^pp 2 ) 

0.745do!o58 ±0.034 

f\l 

i AlVdAi 2 + iyy 2 > 

0.50 ±0.09 ±0.05 

<5|l — S 0 

arg {ApKp) 

1.84 ±0.20 ±0.14 

T%7r7r)o 

IAd7T7r) 0 | 2 /(|A°p 2 + |Ap P | 2 ± |App 2 ) 

0.30 tom ±0.08 

F pf{ 980) 

A p /(98o)| 2 /(|App| 2 + |App 2 + |Ap P 2 ) 

0.29 dp 09 ±0.08 

-^(7r7r)o (7T7r)o 
^-L 7T7r)o 

^_L — bpf( 980) 

^(7T7r)o(7r7r)o ^0 

\A^) 0 (itk) 0 \ 2 /(\A G PP \ 2 ± |A|lp| 2 + App 2 ) 
arg {A^A^J 

ar g(AppA*y (980) ) 

aig(A( 7r7r ) 0 ( W7r ) 0 App) 

0.21 tom ±0.08 
-1.13 ±0.24 

1.92 ±0.24 ±0.16 
3.14 toll ± 0.39 

F 

1 pu 

(|ApJ 2 + |AjL| 2 + |A^J 2 )/(|App|“ + App 2 + |A pp 2 ) 

0.025do!o22 ± 0.020 

rpu) 

■1L 

Ap W | 2 /(|A° w | 2 + |A|L 2 + |A p J 2 ) 

0.70 1^0 ±0.13 

rpUJf 

■'ll 

- % 

l±l 2 /(l±l 2 + \AjS) 
arg (A° pu] A° p * p ) 

0.97 t°of e ± 0-15 
-2.56 t°olt ±0.22 

<*jf - % 

arg(Ap a ,A°*) 

-0.71 toll ±0-32 


arg(A^A; Mo ) 

-1.72 ±2.62 ±0.80 

zpO 

o 

fipf 2 ~ dp(nn) 0 

KJ7(K/ + l±l 2 + iy/) 

ar g( A °/,Kmo) 

0.01 1^2 ±0.03 
-0.56 ± 1.48 ±0.80 

F s+ 

aiiT 

s:S+ JT 

u a 1 TT u p(tt-jt) 0 

Kjya^i 2 +1 Ai 2 + oil 2 ) 

arg «t4*(,, )o ) 

i a +1.0 +1.2 

1 ^ —0.7 -0.8 

-0.09 ++ ±0.38 


the same criteria as those applied to data, following the method described in Ref. 38 


The coefficients are used to determine the efficiency and to build a probability density 
function for each category, which is defined as 


S k (m 1 ,m 2 ,9 l ,9 2 ,<p) = 


Ej= 1 Ei<j 77 e[+A */+*](2 - + X 1 + ViVj)^(m 1 ,m 2 ) 

Ejli J2i<j fte[A^4'](l + ViVj) 


(7) 


The four event categories are used in the simultaneous unbinned maximum likelihood 
fit which depends on the 19 free parameters indicated in Table [3] 

Systematic effects are estimated by fitting with the angular model an ensemble of 1000 
pseudoexperiments generated with the same number of events as observed in data. The 
biases are for the parameters of interest consistent with zero. A systematic uncertainty is 
assigned by taking 50% of the fit bias or the uncertainty on the rms when the latter is 
bigger in order to account for possible statistical fluctuations. 

Several model related uncertainties are envisaged. The B° —s- angular model 

requires knowledge of the lineshape of the af meson. The af natural width is chosen to 
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be 400MeV/c 2 . The difference to the fit results obtained by varying the width from 250 to 
600MeV/c 2 is taken as the corresponding systematic uncertainty. In addition, a systematic 
uncertainty is obtained by introducing the CP-odd component in the fit model of the decay 


amplitude B° —» a^Ti^ by fixing the relative amplitudes of B° —>■ 


a{ 7T 


and P —>• a. 


7T 


components to the values measured in Ref. 39 . Another source of uncertainty originates 


in the modelling of the low mass (7r + 7r _ ) S-wave lineshape. The /o(500) mass and natural 
width uncertainties from Ref. |34| and the uncertainty on the parameter that quantifies the 
non-resonant contribution are propagated to the angular analysis parameters by generating 
and fitting 1000 pseudoexperiments in which these input values are varied according to a 
Gaussian distribution having their uncertainties as widths. The root mean square of the 
distribution of the results is assigned as a systematic uncertainty. The same strategy is 
followed to estimate the systematic uncertainties originating from the p°, / 0 (500) and c o 
lineshape parameters. 

The uncertainty related to the background subtraction method is estimated by varying 
within their uncertainties the fixed parameters of the mass fit model and studying the 
resulting angular distributions and two-body mass spectra. The difference to the fit 
results is taken as a systematic uncertainty. An alternative subtraction of the background 
estimated from the high-mass sideband is performed, yielding compatible results. 

The knowledge of the acceptance model described in Eq.Q comes from a finite sample 
of simulated events. An ensemble of pseudoexperiments is generated by varying the 
acceptance weights according to their covariance matrix. The root mean square of the 
distribution of the results is assigned as a systematic uncertainty. 

The resolution on the helicity angles is evaluated with pseudoexperiments resulting 
in a negligible systematic uncertainty. The systematic uncertainty related to the (7r + 7r _ ) 
mass resolution is estimated with pseudoexperiments by introducing a smearing of the 
(7r + 7r _ ) mass. Differences in the parameters between the fit with and without smearing 
are taken as a systematic uncertainty. 

Table [4] details the contributions to the systematic uncertainty in the measurement of 
the fraction of B° —)■ p°p° signal decays in the B° —> (7r + 7r _ )(7r + 7r _ ) and its longitudinal 
polarisation fraction. 

The final results of the combined two-body mass and angular analysis are shown 
in Figure |3] and Table [ 3 } The fit also allows for the extraction of the fraction of B° —> p°p° 
decays in the B° —> (tt + 7t~)(tt + 71~) sample, defined as 


P(B° ->■ p°p°) 


E}ii Ei<j KelAiAtuij] ’ 


( 8 ) 


which is 

P(B° —> p°p°) = 0.619 ± 0.072 (stat) ± 0.049 (syst). 


The B° —y p°p° signal significance is measured to be 7.1 standard deviations. The 
significance is obtained by dividing the value of the purity by the quadrature of the 
statistical and systematic uncertainties. No evidence for the B° —y p°/ 0 (980) decay mode 
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Table 4: Relative systematic uncertainties on the longitudinal polarisation parameter, /l, and 
the fraction of B° —> p°p° decays in the B° —> sample. The model uncertainty 

includes the three uncertainties below. 


Systematic effect 

Uncertainty on / L (%) 

Uncertainty on P(B° —> p°p°) (%) 

Fit bias 

0.1 

0.8 

Model 

3.6 

6.2 

B° -> ai(1260)+vr- 

1.2 

1.1 

S-wave lineshape 

3.4 

6.1 

Lineshapes 

<0.1 

0.1 

Background subtraction 

0.1 

0.5 

Acceptance integrals 

2.7 

4.5 

Angular/Mass resolution 

0.8 

1.5 



Figure 3: Background-subtracted p, cos$i ,2 and <p distributions. The black dots 

correspond to the four-body background-subtracted data and the black line is the projection of 
the fit model. The specific decays B° — > p°p° (brown), B° — > ujp° (dashed brown), B° — »• VS 
(dashed blue), B° — > SS (long dashed green), B° —> VT (orange) and B° — > afn T (light blue) 
are also displayed. The B° —> p°p° contribution is split into longitudinal (dashed red) and 
transverse (dotted red) components. Interference contributions are only plotted for the total 
(black) model. The efficiency for longitudinally polarized B° —> p°p° events is ~5 times smaller 
than for the transverse component. 


is obtained. The fraction of longitudinal polarisation of the B° -P- p°p° decay is measured 
to be 


/ L = 0.7451^058 ( stat ) ± °- 034 (syst). 
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5 Branching fraction determination 

The branching fraction of the decay mode B° —)■ p°p° relative to the decay B° —y (j)K *(892)° 
can be expressed as 


B(B°-> p°p 0 ^ 


A/ l ■ P(B °—* p°p°) N'(B 0 —)■ (7T + 7T )(7T + 7T )) 


B{B(j)K* (892)°) P(R 0 ^ 0A'*(892)°) 


x 


N'(B 0 ^ (K+K~)(K+ 7T-)) 

P(0 -)■ -)■ A' +7T-) 


X 


P(p° —> 7T+7T - ) 2 


(9) 


where the factor A/ L corrects for differences in detection efficiencies between experimental 
and simulated data due to the polarisation hypothesis of the B° —> p°p° sample, P(B° —* 
p°p°) and P(B °—» 0A'*(892)°) are the fractions of P° —» p°p° and B° —> (j)K*{ 892)° signals 
in the samples of P° —>■ (7r + 7r _ )(7r + 7r _ ) and P° —» (A' + A' _ )(A' + 7 t _ ) decays, respectively. 
The quantities N'(B°—> and N'(B° —>• (A' + A'~)(A" + 7t _ )) are the yields of 

P°—> (7r + 7r _ )(7r + 7r _ ) and P° —> (K + K~)(K + n~) decays as determined from a fit to the 
four-body mass distributions, weighted for each data-taking period by the efficiencies of the 
signal and normalisation channels obtained from their respective simulated data. Finally, 
B{(j> —* A' + A' - ), P(A'*(892)° —> K + tt~ ) and £>(p° —> 7r + 7r ~) denote known branching 
fractions |2l]. 

The product A f L ■ P(B° —> p°p°) is determined from the amplitude analysis to be 
1.13 ± 0.19 (stat) ± 0.10 (syst). This quantity is mainly related to the modelling of the 
S-wave component, and dominates the systematic uncertainty of the parameters of interest. 

The fraction of B° —y (f>K* (892)° present in the B° —» (K + K~)(K + n~) sample is taken 
from Ref. 40 . A 1% systematic uncertainty is added, accounting for differences in the 


selection acceptance for P- and S-wave contributions. 

The amounts of B° —>■ (7r + 7r - )(7r + 7r - ) and B° —>■ (K + K~)(K + tt~) candidates are 
determined from the four-body mass spectra analysis and their associated statistical 
and systematical uncertainties are propagated quadratically to the branching fraction 
uncertainty estimate. 

The limited size of the simulated events samples that meet all selection criteria result 
in a systematic uncertainty of 1.7% (2.6%) on the measurement of the relative branching 
fraction for the 2011 (2012) data-taking period. The impact of the discrepancies between 
experimental and simulated data related to the B° meson kinematical properties is 0.6% 
(1.2%). The efficiencies of the particle-identification requirements are determined from 
control samples of data with a systematic uncertainty of 0.5%, mostly originating from the 
limited size of the calibration samples. An additional 1% systematic uncertainty on the 
tracking efficiency is added accounting for different interaction lenghts between n and K. 

The relative branching fraction is measured to be 


B(B°^p°p 0 ) 
B(B 0 ^ 0AT*(892)°) 


= 0.094 ± 0.017 (stat) ± 0.009 (syst). 


( 10 ) 
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The agreement between the results obtained in the two data-taking periods is tested with 
the best linear estimator technique 41 yielding compatible results. 

The average branching fraction of B° —» 4>K*( 892)° as determined in Ref. 1211 does 
not take into account the correlations between systematic uncertainties due to the S-wave 
modelling. Instead, we average the results from Refs. [42 -44 including these correlations 
to obtain B(B° —* 0A'*(892)°) = (1.00 ± 0.04 ± 0.05) x 10~ 5 . Using this value in Eq.(10), 
the branching fraction of B° —>- p°p° is 


B(B° p°p°) = (0-94 ± 0.17 (stat) ± 0.09 (syst) ± 0.06 (BF)) x 10 


-6 


where the last uncertainty is due to the normalisation channel branching fraction. Using 
the B° —> p°p° branching fraction, the p°/o(980) amplitude, a phase space correction and 
assuming 100% correlated uncertainties, an upper limit for the B° —> p°/ 0 (980) decay, at 
90% confidence level, is obtained 

B{B° ->• p°/ 0 (980)) x £>(/ 0 (980) ->• n + n~) < 0.81 x 10“ 6 . 


6 Conclusions 

The full data set collected by the LHCb experiment in 2011 and 2012, corresponding to an 
integrated luminosity of 3.0fb 1 , is analysed to search for the B°—t p°p° decay. A yield of 
634 ± 28 ± 8 R°—>■ signal decays with 7r + 7r - pairs in the 300-1100 MeV/c 2 

mass range is obtained. An amplitude analysis is conducted to determine the contribution 
from B° —> p°p° decays. This decay mode is observed for the first time with a significance 
of 7.1 standard deviations. In the same 7r + 7r - pairs mass range, B ° —> (7r + 7r - )(7r + 7r - ) 
decays are also observed with a statistical significance of more than 10 standard deviations. 

The longitudinal polarisation fraction of the B° —> p°p° decay is measured to be 
/l = 0.7451q;o58 (stat) ± 0.034 (syst). The measurement of the B° —>■ p°p° branching frac¬ 
tion reads 

B{B° p°p°) = (0-94 ± 0.17 (stat) ± 0.09 (syst) ± 0.06 (BF)) x 10~ 6 , 

where the last uncertainty is due to the normalisation channel. These results are the most 
precise to date and will improve the precision of the determination of the CKM angle a. 

The measured longitudinal polarisation fraction is consistent with the measured 
value from BaBar [8| while it differs by 2.3 standard deviations from the value obtained 
by Belle [9|. The branching fraction measurement is in agreement with the values measured 
by both BaBar |8ij and Belle |[9] collaborations. 

The evidence of the B° —> p°/ 0 (980) decay mode reported by the Belle collaboration |9 
is not confirmed, and an upper limit at 90% confidence level is established 

B(B° ->• p°/ 0 (980)) x £>(/ o (980) Tr+Tr”) < 0.81 x 10" 6 . 
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